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We establish a connection between the cascading of gluon momenta modeled within the diffusive 
approximation of the Balitsky-Fadin-Kuraev-Lipatov kernel and the thermodynamics of directed 
polymers on a tree with disorder. Using known results on the low-temperature spin-glass phase of 
this statistical mechanic problem we describe the dynamical phase space of gluon transverse mo- 
menta near saturation including its fluctuation pattern. It exhibits a nontrivial clustering structure, 
analoguous to "hot spots" , whose distributions are derived and possess universal features in common 
with other spin-glass systems. 



1. Saturation in QCD is expected to occur when parton densities inside an hadronic target are so high that their 
wave-functions overlap. This is expected from the high-energy evolution of deep-inelastic scattering governed by the 
Balitsky Fadin Kuraev Lipatov (BFKL) kernel Q ■ The BFKL evolution equation is such that the number of gluons 
of fixed size increases exponentially and would lead, if not modified, to an average violation of unitarity. Saturation 
\& ■ is characterized by a typical transverse momentum scale Q S (Y), depending on the overal rapidity of the reaction, 

at which the unitarity bound is reached by the BFKL evolution of the amplitude. The problem we want to address 
here is the characterization of the gluon-momentum distribution near saturation. This study may help analyzing the 
Color Glass Condensate (CGC) which has been proposed 0] as a description of the QCD phase at full saturation. As 
we shall see, our aim in the present paper is to fully understand the spectrum of transverse momentum fluctuations 
| among the gluons which are generated by the BFKL evolution in rapidity, near saturation. 

Indeed, one of the major recent challenges in QCD saturation is the problem of taking into account the role of 
I/"") \ fluctuations which are present in the BFKL evolution, i.e. in the dilute regime. Quite surprisingly, the fluctuations in 
this region modify the overall approach to saturation by a mechanism which has been understood as the formation of 
traveling-wave solutions of nonlinear differential equations. If at first one neglects the fluctuations (in the mean-field 
approximation), the effect of saturation on a dipole-target amplitude is described by the nonlinear Balitsky-Kovchegov 
! [El (BK) equation, where a nonlinear damping term adds to the BFKL equation. As shown in 6] , this equation falls 
into the universality class of the Fisher and Kolmogorov Petrovsky Piscounov (F-KPP) nonlinear equation which 
admits asymptotic traveling- wave solutions of "pulled-front" type Q ■ In the QCD case, this means that they are driven 
by the gluons of higher momentum. Hence the overall behaviour of the solution is driven by the high-momentum tail, 
which is precisely in the "dilute" domain of the BFKL regime, far from the range which is dominated by saturation. 
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The exponential behaviour of the BFKL evolution quickly enhances the effects of the tail towards a region where 
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finally the nonlinear damping regulates both the traveling-wave propagation and structure. 

In these conditions it was realized, for "pulled-front" traveling waves [t| and in QCD [lCj, that the fluctuations 
inherent of the dilute regime may have a surprisingly large effect on the overall solution of the nonlinear equations 
of saturation. Indeed, a fluctuation in the dilute regime may grow exponentially and thus modify its contribution to 
the overall amplitude. Hence, in order to enlarge our understanding of the QCD evolution with rapidity, it seems 
• • ■ important to give a quantitative description of the fluctuation pattern generated by the BFKL evolution equations 
. ' for the set of cascading dipoles (and thus gluons) near saturation, which is the goal of our paper. We shall work in 
■ the leading order in l/N c , where the QCD dipole framework is valid [ill ll^. 

The equations for the generating functionale of moments of the QCD dipole distributions have been derived already 
some time ago |l lj and they have been used in instructive numerical simulations of onium-onium scattering [Llj ] and 
also applied recently in diffractive and non-diffractive deep- inelastic scattering |l4j . However, to our knowledge, a 
thorough investigation of their solutions for the pattern of dipoles or equivalently for the event-by-event properties of 
gluon transverse momenta has not yet been considered. We will provide an answer to this question in the QCD dipole 
approximation (and moreover in the diffusive approximation of the 1-dimensional BFKL kernel, see further on), and 
the results appear quite rich already within this approximation scheme. Many aspects we find show "universality" 
features and thus are expected to be valid also in the dipole framework beyond the approximations (going beyond the 
dipole approximation is still a conceptual challenge [la])- 

The strategy of the present paper is to determine the structure of the event-by-event gluon transverse-momentum 
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distributions. In the present paper, and as a first step towards a more complete solution, we will simplify the problem, 
focusing on modeling the QCD gluon cascading using the diffusive approximation of the BFKL kernel and aiming at 
a quantitative understanding of the transition-to-saturation phase within this scheme. 

We show that one can build a mapping between the energy evolution of the set of gluon momenta (in the 1- 
dimensional approximation) and the time evolution of directed polymers on a tree with disorder. This is obtained 
through the matching of the partition function of the polymer problem with a related event-by-event physical quantity 
defined from the cascading gluon process. Then we use the relation between the cascading evolution around the 
unitarity limit of the amplitude and the saturation mechanism to find that a spin glass structure of the gluon phase 
space exists at saturation, in direct analogy with the low-temperature phase of the directed polymers. We then derive 
the physical consequences of the spin-glass phase for the distribution of gluon transverse momenta at saturation in 
our model of cascading. 

The plan of our study is the following. In section 2, we shall recall the main ingredients of the rapidity evolution of 
the cascading dipoles and gluon-momentum distributions as driven by the BFKL evolution in the 1-dimensional and 
diffusive approximations. In section 3, we introduce the mapping between the cascading of gluons and the propagation 
of directed polymers on random trees. In section 4 we recall the nontrivial properties of the directed polymer problem, 
characterized by a spin- glass transition at low temperature. In section 5, we provide a check of the consistency of the 
gluon/polymer mapping with the saturation mechanism in the mean field approximation. In section 6, we derive the 
gluon phase space properties at large rapidity from the spin-glass features. The final section 7 is devoted to conclusion 
and outlook. 

2. In this section, we present some known results on the equations governing the dipole distributions in the BFKL 
framework together with those concerning the saturation amplitude in the mean-field approximation and relevant for 
our study. 



QCD rapidity evolution 

Let us start by briefly describing the QCD evolution of the dipole distributions in the Balitsky-Fadin-Kuraev-Lipatov 
BFKL framework [H 111 ■ 

Formally, let us denote Z vw [Y, U] the generating functional of moments for the N-uple dipole-probability distri- 
butions after a rapidity evolution range Y. The dipole probabilities can be deduced from Z vw [Y, U] by functionally 
differentiating with respect to the source fields U and then letting U — > 0. the evolution equation satisfied by Z VW [Y, U] 
reads 
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is the BFKL kernel 0] describing the dissociation vertex of one dipole (v,w) into two dipoles at (v,z) and (z,w), 
where v, w, z are arbitrary 2-dimensional transverse space coordinates . 

From Eq.Q, we will mainly retain its physical meaning. The functionale structure of Eq.JQ) describes a 2- 
dimensional tree structure of dipoles in transverse position space evolving with rapidity. Let us, for instance, focus 
on the rapidity evolution starting from one massive qq pair or onium |lll| see Fig^ At each branching vertex, the 
wave function of the onium-projectile is described by a collection of color dipoles. The dipoles split with a probability 
per unit of rapidity defined by the BFKL kernel J2J- The BFKL rapidity evolution generates a cascade of dipoles. 
Each vertex of the cascade is governed in coordinate space by the 2-dimensional BFKL kernel J2J), determining both 
the energy and transverse position space evolution of the cascade, as sketched in FigQ] As we shall see later on, this 
cascading property can be converted in a similar cascade of gluons in momentum space. 

At this stage, it is interesting to note the mathematical similarity of Q with the BK equation as verified 
by the S-Matrix element of the dipole-target amplitude S(v,w,Y) in the mean-field approximation. The equation is 
indeed formally the same if one replaces term by term the functional Z VW [Y~, U] by the function <S(v, w, Y) representing 
the S-Matrix element. 
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FIG. 1: BFKL cascading and saturation. The QCD branching process in the BFKL regime and beyond is represented along 
the rapidity axis (upper part). Its 2-dimensional counterpart in transverse position space is displayed at two different rapidities 
(lower part). The interaction region is represented by a shaded disk of size 1/Q. At rapidity Y\, the interaction still probes 
individual dipoles (or gluons), which corresponds to the exponential BFKL regime. There exists a smooth transition to a regime 
where the interaction only probes groups of dipoles or gluons, e.g. at rapidity Y2. This gives a description of the near- saturation 
region corresponding to a mean- field approximation, where correlations can be neglected. Further in rapidity, Y > Y2, other 
dynamical effects, such that merging and correlations appear, leading eventually to the CGC. 



Gluons in the 1-dimensional and diffusive approximation 

As an approximation of the 2-dimensional formulation of the BFKL kernel obtained when one neglects the 
impact-parameter dependence, we shall restrict our analysis in the present paper to the 1-dimensional reduction of 
the problem to the transverse-momenta squared k\ of the cascading gluons. After Fourier transforming to transverse- 
momentum space, the leading-order BFKL kernel Q defining the rapidity evolution in the 1-dimensional approxima- 
tion is known |a] to act in transverse momentum space as a differential operator of infinite order 

x (-9 i ) = 2^(l)-V(-^)-V(l + ^) (3) 

where I = log k 2 and Y is the rapidity in units of the fixed coupling constant as ^ a . 

In the sequel, we shall restrict further our analysis to the diffusive approximation of the BFKL kernel. We thus 
expand the BFKL kernel to second order around some value j c 

Xil) ~ Xc + x'cil ~ 7c) + IXc(7 " 7c) 2 = 4) - 47 + Ml 2 , (4) 
where 7 C will be defined in such a way to be relevant for the near- saturation region of the BFKL regime. 



Within this diffusive approximation, it is easy to realize that the first term (Aq) is responsible for the exponential 
increase of the BFKL regime while the third term (A 2 ) is a typical diffusion term. The second term (A% ) is a "shift" 
term since it amounts to a rapidity-dependent redefinition of the kinematic variables, as we shall see. 

In Eq.(@J, 7 C is chosen in order to ensure the validity of the kernel Q in the transition region from the BFKL 
regime towards saturation. Indeed, the derivation of asymptotic solutions of the BK equation [6j leads to consider 
the condition 

X(lc) = 7c x'ilc) (5) 

whose solution determines j c . 

This condition applied to the kernel formula © gives j c = y/ Aq /A2 ~ 0.6275... and {Ao, Ai, A2} « 
{9.55,25.56,24.26}. It is worth noting that a more general choice of the Aj's can be considered for the forthcom- 
ing discussion, leaving the possibility of considering also "effective" kernels eventually taking into account next 
leading corrections and thus other values of j c . 



Saturation in the mean-field approximation 

As already noticed 14], the 2-dimensional BK equation describing saturation in the mean-field approximation is 
formally similar to Q by substituting Z by S, the S-matrix element of the scattering amplitude. In the 1-dimensional 
approximation which we will use, the BK equation reads 



dyT(l, Y) = x(-dt) T(l, Y) - T 2 (l, Y) 



(6) 



with the same notations as in @. 

In the 1-dimensional approximation, the Fourier transform of the dipole-target scattering amplitude T(l,Y) is 
related to the unintegrated gluon distribution in the target (by a simple convolution). Applying only the linear 
operator J2J to T(l, Y) leads to an exponential increase of the gluon density with rapidity characteristic of the BFKL 
regime. The damping term in Eq.^HJ is responsible for the saturation mechanism in the mean-field approximation. 

There exists an already long literature on the properties of the BK equation (see, for instance, 0, Hila] ) • For further 
purpose, we remark that many results which can be obtained from the solutions of the BK equation can already 
be derived from the BFKL evolution when only the linear operator in Eq.JBJ) is retained. The idea [lil is to 
consider the extension of the BFKL evolution in the rapidity region when the unitarity limit is reached and assume 
a smooth transition to saturation. This amounts to study the BFKL regime alongside the saturation line defined by 
the saturation scale Q S (Y) in the 2-dimensional k 2 ,Y plane. 

Indeed, one obtains a significant part |l8l [19( of the results obtained |fj by a direct derivation of the asymptotic 
solutions of the BK equation. Let us list some of the most relevant results of this approach based on the unitarity 
bound: The gluon momentum distribution can be derived at high rapidity and large enough k. One finds 
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The saturation scale Q 2 (Y) is given by 
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Furthermore, at large rapidity, the third term of Q and the subasymptotic correction in (JHJ) become negligible, leaving 
a gluon momentum distribution exhibiting geometric scaling [20| 
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The asymptotic results (j^Jl can be directly obtained from the unitarity bound applied to the solution of BFKL equation 
[Isj . while the subasymptotic results can also be obtained from the BFKL equation, but with more care taken about 
the boundary conditions [l9j . 



3. Now we introduce the mapping from the BFKL regime of QCD to the thermodynamics of directed polymers. 
As we have mentionned earlier, the BFKL rapidity evolution generates a cascade of gluons ^lj- Each vertex of the 



cascade is governed in coordinate space by the 2-dimensional BFKL kernel, determining together the energy and the 
transverse space evolution of the cascade, see Fig^ In transverse-momentum space and within the 1-dimensional 
diffusive approximation we already noticed that the BFKL kernel models boils down to a branching, shift and 
diffusion operator acting in the gluon transverse-momentum-squared space. Hence the cascade of gluons can be put 
in correspondence with a continuous branching, velocity-shift and diffusion probabilistic process whose probability by 
unit of rapidity is defined by the coefficients Ai of (@J. 

Let us first introduce the notion of gluon- momenta "histories" ki{y). They register the evolution of the gluon- 
momenta starting from the unique initial gluon momentum fc(0) and terminating with the specific i-th momentum 
hi, after successive branchings. They define a random function of the running rapidity y, with < y < Y, the final 
rapidity range when the evolution ends up (say, for a given total energy). It is obvious that two different histories 
ki(y) and kj(y) are equal before the rapidity when they branch away from their common ancestor, see Fig^ 

We then introduce random paths Xi(t) using formal space and time coordinates x, t , < t < At which we relate 
to gluon-momenta "histories" as follows: 

y=4~; logkUy) = -/3(xi(t)-x(0)) + (A -A 1 )y (10) 

^0 

where (Aq — A\)y is a conveniently chosen and deterministic "drift term", a;(0) is an arbitrarily fixed origin of an 
unique initial gluon and thus the same for all susequent random paths. The random paths Xi{t) are generated by a 
continuous branching and Brownian diffusion process in space-time (cf. FigEJ - 

As we shall determine later on, the important parameter /?, which plays the role of an inverse temperature for the 
Brownian process, will be fixed to 



In fact, the relation (|ll|l will be required by the condition that the stochastic process of random paths describes 
the BFKL regime of QCD near saturation. Another choice of /3 would eventually describe the same branching process 
but in other conditions. Hence the condition (| 1 1 |t will be crucial to determine the QCD phase at saturation (within 
the diffusive approximation). 

Let now introduce the tree-by-tree random function defined as the partition function of the random paths system 

Z{t) = J2 e~ Px ' W = e-^ ^ x ^= l ^ V) • (12) 

Using the fact that the total number of momenta in one event at rapidity y is n w ( n ) = e A " y at large y, (I12|) takes 
the form 

Z = A e Aiy x k 2 (y) , (13) 

where ~ Y27=i ^i(v) = k 2 {y) is the event-by-event average over gluon momenta at rapidity y and A = e~@ x ° is an 
arbitrary constant. Note that one has to distinguish ( ••• ), i.e. the sample-to-sample average, from which 
denotes the average made over only one event. 

Z(t) is an event-by-event random function. The physical properties are obtained by averaging various observables 
over the events. Note that the distinction between averaging over one event and the sample-to-sample averaging 
appears naturally in the statistical physics problem in terms of "quenched" disorder: the time scale associated with 
the averaging over one random tree structure is different (much shorter) than the one corresponding to the averaging 
over random trees. 

In parallel with the statistical physics problem |2l| . it is convenient to introduce the generating functional G(t,x) 
of its event-averaged moments ( Z p ) , namely 

G(t,x) = J2 ^ ( [-e ^ Z(t)] P ) = ( exp [-e~^ Z(t)} ) = < exp [-A P(y)/ K 2 } ) , (14) 

where 

y = 4~, log^ 2 = /3x + 4 i ^ ■ (15) 



The moments ( Z p (t)) are easily recovered by differentiation over x or logK 2 at fixed t. 

It is interesting to note the final form of (|14|l showing a similarity (l^j with a unitary S-Matrix amplitude if k is 
interpreted as the scale of a transverse-momentum probe, such as depicted in Fig^ This similarity will be useful for 
the derivation of saturation properties from the BFKL regime in the unitarity limit. 

4. Now that we have established the connection (|10|l between the model of gluon-momenta submitted to branching, 
shift and diffusion, and the random paths defined by #i(t), let us rederive the properties of QCD gluon cascading near 
saturation from the point-of-view of the statistical random-path system. For this sake, we recall the corresponding 
model for directed polymers on a tree, as defined in R.ef . |2lj . 

Directed polymers on a random tree are modeled, see Fig|2 by the time evolution of random paths in 1-dimensional 
space with continuous diffusion and branching (with probability per unit of time fixed respectively to J and 1, by 
convention) . 

Indeed, contemplating the definitions 1|12[) and ljl4Jl it is straightforward to identify Z as the sample-to-sample 
partition function and G(t, x) as the generating fonctional of its sample-averaged moments for the model [2lj | of 
directed polymers on a random tree. 
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FIG. 2: Branching diffusion model for polymers. The coordinates xi(t) ■ ■ ■ xj{t) correspond to the random paths along the tree 
in the x,t phase-space. The oblique axis is for L = /3x + Ai/Ao t, which takes into account the "time-drift" in the mapping to 
the QCD problem. 



Let us list some relevant properties of the polymer problem |21| . 

• u(t, x) = 1 — G(t, x) verifies the Fisher and Kolmogorov Petrovsky Piscounov (F-KPP) equation 

d . 1 d 2 u(t,x) o 
— u(t, x) = - — — h u(t, x)-u (t, x) . 

u(t, x) admits traveling-wave solutions at large time t 0| 

u(t, x) ~ u(x — mp(t)) sa u(x — c(f3)t) , 



(16) 



(17) 



where the wave speed c(/3) is governed by the initial condition u(to,x) ~ e~^ x at large x. The range of values 
of c(/3) = 1/(5 + (3/2 is bounded from below by c{(3 c ) = c = y/2, where the wave velocity remains fixed whenever 
(3 > (3 C = c and (3 C are usually called, respectively, the critical velocity and the critical slope of the waves. 
In the "supercritical" regime (3 > j3 c , the wave front and its speed remain "frozen" and 



m/j(t) = c(*--Iogt) + 0(l) 



• The free energy of the polymer system at equilibrium, defined by 
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has a "frozen" density Fjt — > — c for /3 > [3 C — \/2 indicating the existence of a phase transition at (3 C . This is 
related to the dynamical transition c(/3) — > c in the traveling- wave properties 

• Introducing the tree-by-tree free energy spectrum 

/(i) = -~logi?(t) + ~<logZ(i)) , (20) 

its probability distribution Vp(f) around the average free energy i* 1 is found by deriving asymptotic expressions 
for the moments |2l| to satisfy 

W) ~ -fe^ for /--oo 

W) ~ e^ 2 "^ for / - oo , (21) 

where we quoted only the distribution obtained for the supercritical regime > (3 C . This means that the 
fluctuations of the free energy around its event average remain finite even at asymptotic time, i.e. there is no 
thermodynamical self-averaging, demonstrating that the phase transition is not similar to an ordinary liquid-gas 
phase transition. 

• The corresponding phase transition is towards a spin-glass phase characterized by a non-trivial structure in 
"valleys" or ground state minima, which can be explored quite in detail. We shall describe this structure in 
the next section, since it will give, through matching with the gluon cascading problem, a description of the 
gluon-momentum phase at saturation within the diffusive mean-field scheme. 

5. Let us make now the connection of the directed-polymer properties with the description of the gluon-momentum 
phase near saturation. 

Our aim is to study the structure of the gluon-momentum spectrum when the system reaches saturation. We 
thus start with a simple initial state (e.g. an onium or one gluon) and, the rapidity increasing, the state develops 
into a cascading structure of gluons governed by the model that we assume. It consists in random tree structures 
generated by branching, shifting and diffusing, as previously described. As we have discussed in introduction, the 
rapid exponential increase of the number of gluons generates reinteractions which damp the exponential increase and 
cause the saturation phenomenon Q . 

In the simple mean- field approximation which neglects the dynamical correlations ||] , the saturation mechanism is 
known to be studied from the BFKL properties in the following approximation |l8j . When the exponential increase of 
the gluon density reaches the unitarity limit of the amplitude, the saturation effects set in. This defines a kinematical 
region in the k,Y plane around the saturation scale Q S (Y) where the BFKL regime is expected to smoothly connect 
the saturation region. Indeed, this expectation has been verified for the amplitude, where many aspects of the 
saturation features rigorously found |6j for the solutions of the mean field BK saturation equation have been already 
derived from the BFKL regime considered near the unitarity limit 18] . 

We shall now try and extend this approach of saturation, which provided already valuable results on averaged 
observables, to the event-by-event spectrum of gluon transverse momenta. In other words, we study momenta of the 
order of the saturation scale, i.e. kf = 0(Q 2 S ). This is nothing else than the neighbourhood of the kinematic region 
where "universal" traveling- wave symptotic solutions of the BK equation exist |y] . We thus expect from the studies 
of Ref . that the properties we shall derive will be independent of initial conditions (provided they are leading to a 
"supercritical" regime) , of the precise form of the kernel (beyond the diffusive approximation) and even of the precise 
form of the non-linearities || . 

A comment of caution is nevertheless in order here, concerning the domain of validity of our appoach. The properties 
in the region of smaller momenta ki/Q s 1 will probably depend more on the specific structure of the BFKL kernel 
beyond the diffusive approximation while the region with ki/Q s 3> 1 will keep the memory of the initial conditions. 
They should require a different treatment, which is not included in our present work and deserves later study. 

Inserting the QCD definitions (|13|l in (|19|l for observables after a final rapidity evolution-range Y (or At in time) 
gives 



log/c 2 ) = ( logZ ) — logA — A{Y — (3 mp(At) — A{Y — log A 
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A Y - 7 logF I - A X Y + 0(1) ~ logQ'i ■ (22) 



The second line of Eas. (|22|l is obtained by identifying the inverse temperature j3 with its value defined by the relation 
This the condition to verify k 2 w Q 2 , the identification coming from formula (JSJ. This condition is crucial to 
check the main saturation property that the sample averaging gluon momentum stays at the value of the saturation 
scale. Hence, the relations (122(1 . important for the QCD saturation problem, determines the equivalent temperature 
of the polymer system. Starting from (0J and the inverse temperature obtains /3 = ~ 2.25 > (3 C = y/2 which 
is larger than the critical value. As a consequence, we are led to consider and study the low-temperature phase of the 
equivalent random-path system in order to get information on the QCD phase of gluons near saturation. 

As a further check of the validity of the random cascading model in the "unitarity limit" region, where the average 
momentum is consistently found to be determined by the saturation scale, let us discuss the distribution around the 
average. 

The event- by-event mean momentum k at rapidity Y can be read off from Ea. l|2U[l as being 



= -0f(AD = -^-^f(A o Y) (23) 
and, from the first line of 121(1 . when ~ k 2 /Q 2 (Y) is large, its probability distribution is given by 
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(24) 



By comparison with formula (0, the interpretation of Eq. 1(24(1 is clear. The probability law 1(24(1 for the event-by- 
event partition function coincides with the gluon distribution in the saturation region and follows geometric scaling. 
This gives the confirmation, through the properties of the directed polymers, that the gluon-cascading model near the 
unitarity limit gives the correct saturation scale and gluon distribution. In mathematical terms, this comes from the 
fact that the equation for the generating function of moments for the model ((16(1 is in the same (FKPP) universality 
class than the BK equation 

Both properties 1(22 (1 . 1(24 (1 prove the consistency of the model with the properties expected from saturation. We 
shall then look for other properties of the cascading gluon model. It is important to realize that this consistency 
fails for any other choice of the parameter (3 different from 1(11(1 . This justifies a-posteriori the identification of the 
equivalent temperature of the system in the gluon/polymer mapping framework. 

6. Let us now come to the main topics concerning the determination of structure of the gluon-momentum phase at 
saturation in the diffusive model. 

The striking property of the directed polymer problem on a random tree is the spin-glass structure of the low 
temperature phase. As we shall see this will translate directly into a specific clustering structure of gluon transverse 
momenta in their phase near the "unitarity limit". Let us divide the discussion in a few subsections. 

i) The equivalent temperature. Following the relation ((11(1 , one is led to consider the polymer system at temperature 
T with 



T c (3 
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where j c is the critical exponent. We are thus naturally led to consider the low-temperature phase (T < T c ), at some 
distance T c — T from the critical temperature T c =i/V2. In the language of traveling waves this corresponds to a 
"pulled- front" condition with "frozen" and "universal" velocity and front profile. 

ii) The asymptotic phase landscape. As derived in [21|, the phase space of the polymer problem is structured in 
"valleys" which are in the same universality class as those of the Random Energy Model (REM) 22] and of the infinite 
range Sherrington-Kirkpatrick (SK) model [23| . 

Translating these results in terms of gluon momenta (in modulus), the phase space landscape consists in event- 
by-event distribution of clusters of momenta around some values fc^ = l/(rii) ^2i^ S ikf, where m is the cluster 
multiplicity. The probability weights to find a cluster Sj after the whole evolution range Y is defined by 



W» = 



■ kf 



J2 t k 2 



where the summation in the numerator is over the momenta of gluons within the sf 1 cluster (see FigEJ) - The normalized 
distribution of weights W S i thus allows one to study the probability distribution of clusters. The clustering tree 
structure, (called "ultrametric" in statistical mechanics) is the most prominent feature of spin-glass systems |24j . 



Note again that, for the QCD problem, this property is proved for momenta in the region of the "unitarity limit", or 
more concretely in the momentum region around the saturation scale. This means that the cluster average-momentum 
is also such that fc^ = 0(Q 2 S ). Hence the clustering structure is expected to appear in the range which belongs to the 
traveling-wave front [6j or, equivalently, of clustering with finite fluctuations around the saturation scale. 
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FIG. 3: The clustering structure. The drawing represents the si • • • S4 clusters near momenta k s \ ■ ■ ■ k 3 4. They branch either 
near t <C 1 or (At -i) C 1, where At is the total amount of time evolution. 

iii) The "overlap" function. There exists in statistical physics of spin-glasses a well-known function which reveals 
the space-time structure of a spin-glass state in a quantitative way: the "overlap" function |24| . Translating the 
definitions (cf. |2l| ) in terms of the QCD problem, one introduces the overlap between two final gluon momenta 
ki and kj, using their "histories" ki(y) and kj(y). We recall that ki(y) and kj(y) defined as function of the running 
rapidity y correspond to the two (not always disjoint) paths histories which finally lead to the momenta ki and kj 
observed at final (and large) rapidity Y. 

The overlap is defined by the quantity Qij, the fraction of the total rapidity Y for which ki(y) — kj(y), < Qij < 1 . 
It is a measure of the rapidity interval ( "time" ) during which they belong to the same branch before splitting. For a 
given event, the probability y(q) of finding an overlap q can be defined by 

1 ™ 

y(q) dq = V kj k] Q(q < Qij < q + dq ^ h{qY) = k {qY) in [q,q + dq}) , (26) 

{L(X)} i,i=i 

where the characteristic function O restricts the sum to thepairs of momenta which where on the same branch in the 
interval of rapidity [q,q + dq]. The result reads (from Ref.|2j|) 

y(q) dq =S(q-l)y + S(q) (1 - y) = S(q 1) £ Wl + S(q) ( 1 - ]T WU , (27) 

si \ si / 

where y = J2si ^si ^ s an event-by-event indicator of the strength of clustering. 

The result H27J) means that, in the thermodynamical limit, the gluon momenta branch effectively either near y = 
or near y — Y for asymptotic Y , see Figdl Hence the "microscopic" branching structure of the model of FigElat finite 
y, occupies only a vanishing (when Y — > 00, i.e. in the thermodynamical limit) fraction of the rapidity evolution, 
either at the beginning or at the end but not in the middle nor for a finite fraction of the asymptotic rapidity. 
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FIG. 4: The probability distribution of overlaps FL(y). The figure (the simulation is by courtesy from j^J, using the method 
of Ref. 26]) is drawn both for the theoretical QCD value 7 C = 0.7, and for y c = 0.3 for comparison. The statistics used for the 
simulation is 10 s events in 10 3 bins for y c = 0.3 and 3.25 10 6 events in 10 2 bins for 7 C = 0.7. 



The probability distribution of overlaps 11(3^) is identical to the one of the REM and SK models and shares many 
qualitative similarities with other systems possessing a spin-glass phase 25, 26]. Being quite non-trivial analytically, 
it is more easily defined by its moments 

ir >_ <-£^ f^»— for .-1 <,<„, (28, 

T{2v) T(n-v) Jo dp 

poo 

g(j>)= dvv- T ' T '- 1 (l-e-'-^J . (29) 



with (see 122, 26]) 



The overlap distribution depends only on T/T c , which confirms its high degree of universality . One finds, in particular, 
( y ) = 1 — T/T c , which is the average probability of having overlap 1, while T/T c is the probability of overlap 0. 

The rather involved probability distribution IL(y) is quite intringuing, as obtained by inversion of moments ( y v ) 
|25l |2f|. It possesses apriori an infinite number of singularities at y = l/n,n integer. It can be seen when the 
temperature is significantly lower from the critical value, see for instance the curve for 1 — T/T c = .7 in FigQ] 
However, the predicted curve for the QCD value 1 — T/T c = .3 is smoother and shows only a final cusp at W s = 1 
within the considered statistics. It thus seems that configurations with only one cluster can be more prominent than 
the otherwise smooth generic landscape. However, it is also a "fuzzy" landscape since many clusters of various sizes 
seem to coexist in general. 

iv) The cluster distribution. The sample-to-sample distribution of "valleys" can be explicitely computed by a 
generalization of the overlap function |2fil l26| to any power e, W^, which translates into a well-defined distribution 
of cluster-weight probability. For instance one determines: 



W -~t.-i (i _ w y, 

r( 7c )r(i - 7c ) 



HWuW 2 ) = 7c ( 1 %, 9 \ r2n 1 P 1 x6(l-W 1 -W 2 ) 



(WiW 2 ■ ■ ■ Wp)-*'- 1 (i - Ei ^) P7c 



w 2 , — ,w p ) = 7? r(p) — ;r v^r V 1 '"' — xW-Y, w i) . ( 3 °) 



r(p 7c )rp(i-7 C ) 



for, respectively, the average number $(VF) of clusters of weight W S i = W, the average number of pairs $(Wi, W2) of 
clusters in the same event with weights W S i = Wi, W2, and more generally W2, ■ ■ ■ , W p ) the average number 

of p-uples of clusters with weights W S i = W2, ■ ■ ■ W p . 

It is interesting to note that smaller is -f c — T/T c , smaller is the probability to find W s ~ 0, i.e. no clustering, and 
higher is the probability for having one big cluster W s ~ 1. As an illustration, the phenomenological value 7 C ss 0.3 
leads to a more pronounced clustering structure than the BFKL value 7 C « 0.6. The sharp structuration in clusters 
expected for small j c can already be seen in the overlap distribution of Fig^J 

7. Let us summarize the main points of our approach. We investigated the landscape of transverse momenta in 
gluon cascading around the saturation scale at asymptotic rapidity. Limiting our study to a diffusive 1-dimensional 
modelization of the BFKL regime of gluon cascading, we make use of a mapping on a statistical physics model 
for directed polymers propagating along random tree structures at fixed temperature. We then focus our study 
on the region near the unitarity limit where information can be obtained on saturation, at least in the mean-field 
approximation. Our main result is to find a low-temperature spin-glass structure of phase space, characterized 
by event-by-event clustering of gluon transverse momenta (in modulus) in the vicinity of the rapidity-dependent 
saturation scale. The weight distribution of clusters and the probability of momenta overlap during the rapidity 
evolution are derived. 

Interestingly enough the clusters at asymptotic rapidity are branching either near the beginning ("overlap 0" or 
y/Y <C 1) or near the end ("overlap 1" or 1 — y/Y <C 1) of the cascading event. The probability distribution of 
overlaps is derived and shows a rich singularity structure. 

Phenomenological and theoretical questions naturally emerge from this study. 

On a phenomenological ground, it is remarkable that saturation density effects are not equally spread out on the 
event-by-event set of gluons; our study suggests that there exists random spots of higher density whose distribution 
may possess some universality properties. In fact it is natural to expect this clustering property to be present not only 
in momentum modulus (as we could demonstrate) but also in momentum azimuth- angle. This is reminiscent of the 
"hot spots" which were some time ago advocated from the production of forward jets in deep-inelastic scattering 
at high energy (small Bjorken variable). The observability of the cluster distribution through the properties of "hot 
spots" is an interesting possibility. 

On the theoretical ground, the questions arise how the pattern of fluctuations we found influence the onset of 
saturation and its properties, given the fact that fluctuations are of primordial importance in the development of the 
whole QCD process, as we have stressed in the introduction. 

First, the question arises to extend the study to the full BFKL kernel and in two dimensions, as given initially 
by Eq.Q. Universality properties of the solutions of the BK equations in the traveling- wave formalism |6j and 
in the 2-dimensional case |22( would encourage such a generalization. However, some care must be taken of the 
potential problems which can arise from an event-by-event description of the cascading regime. The suggestion is 
[2jj to formulate the problem directly in the full 2-dimensional position space, for which the gluon cascading and the 
unitarity constraints are well defined. Some preliminary numerical results are encouraging. 

Second, it would be instructive to incorporate the pattern of fluctuations we found in the calculation of the contri- 
bution of Pomeron loops to the scattering amplitude. The Pomeron loops are found to be important to consider and 
are the subject of many studies (see, e.g. f 0J). The Pomeron loops are obtained by merging the splitting of cascading 
gluons of the projectile with the merging in the target. We think that the spin-glass structure of the cascading gluons 
could lead to useful hints for dealing with the Pomeron loops. 

In conclusion, we think that the already rich structure of the phase space landscape of gluon clusters found with 
the diffusive model of gluon cascading opens an interesting way towards a deeper understanding of the QCD phase 
in the saturation regime and could lead to a deeper description the Color Glass Condensate. 

I would like to thank Edmond Iancu for inspiring and stimulating remarks and Cyrille Marquet and Gregory Soyez 
for many fruitful suggestions and a careful reading of the manuscript. Special thanks are due to Gregory Soyez for 
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